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Abstract 

Some properties of the polarized neutron matter at finite temperature has been studied using the 
lowest order constrained variational (LOCV) method with the AVis potential. Our results indicate 
that spontaneous transition to the ferromagnetic phase does not occur. Effective mass, free energy, 
magnetic susceptibility, entropy and the equation of state of the polarized neutron matter at finite 
temperature are also calculated. A comparison is also made between our results and those of other 
many-body techniques. 
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I. INTRODUCTION 



The spontaneous phase transition to a ferromagnetic state in the neutron matter is of par- 
ticular interest in astrophysics. Specially, this transition could have important consequences 
for the physical origin of magnetic field of pulsars, that are believed to be rapidly rotating 
neutron stars with strong surface magnetic fields in the range of 1012-1013 Gauss 111, fly, 3, 
and also the evolution of a protoneutron star. A hot neutron star is born within a short 
time after supernovae explosion. In this stage (protoneutron star), the interior temperature 
of a neutron star is of the order 20-50 MeV [sj. Therefore, the study of magnetic properties 
of polarized neutron matter at finite temperature is of special interest in the description of 
protoneutron stars. 

There exist several possibilities of the generation of the magnetic field in a neutron star. 
From the nuclear physics point of view, one is the possible existence of a phase transition to a 
ferromagnetic state at densities corresponding to the theoretically stable neutron stars and, 
therefore, of a ferromagnetic core in the liquid interior of such compact objects. Such a pos- 
sibility has been studied by several authors using different theoretical approaches [6-29] , but 
the results are still contradictory. In most calculations, neutron star matter is approximated 
by pure neutron matter, as proposed just after the discovery of pulsars. Whereas some 
calculations, like those of based on hard sphere gas model Q], Skyrmelike interactions 



23l |. Reid soft-core potential [13| and relativistic Dirac-Hartree-Fock approximation with 



an effective nucleon-meson Lagrangian 181] showed that neutron matter becomes ferromag 



221 ] and Brueckner-Hartree-Fock 



netic at some densities. Others, like recent Monte Carlo 
calculations [21-23] using modern two-body and three-body realistic interactions show no 
indication of ferromagnetic transition at any density for neutron matter and asymmetrical 
nuclear matter. 

Most of the studies of the ferromagnetic transition in neutron matter and nuclear matter 
have been done at zero temperature. The properties of polarized neutron matter both at 



zero and finite temperature have been studied by several authors [29|, 3^,1311] • Bombaci et al. 
(BPRRV) 31| have studied the properties of polarized neutron matter within the framework 
of the Brueckner-Hartree-Fock formalism using the AVis nucleon-nucleon interaction. Their 
results show no indication of a ferromagnetic transition at any density and temperature. 
Lopez- Val et al. 301] have used the Dl and DIP parameterization of the Gogny interac- 
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tion and the results of their calculation show two different behaviors: whereas the DIP 
force exhibits a ferromagnetic transition at density of p ~ whose onset increases 

with temperature, no sign of such a transition is found for Dl at any density and temper- 



ature. Rios and Polls 



2^ have used Skyrme-like interactions and their results indicate the 



occurrence of a ferromagnetic phase of neutron matter. 



Recently, we have computed the properties of pola^zed ueutron matter M. pol-ized 

symmetrical nuclear matter [33!], such as total energy, magnetic susceptibility, pressure, etc at 
zero temperature using the microscopic calculations employing the lowest order constrained 
variational (LOCV) method with the AVis potential. We have also calculated the properties 
of spin polarized asymmetrical nuclear matter and neutron star matter 3J] using the LOCV 



method employing the AVis 



35|, Reid93 [36], UV^ |37| and AVi^ |38| potentials. We have 



concluded that the spontaneous phase transition to a ferromagnetic state in the neutron 
matter, symmetrical and asymmetrical nuclear matter and neutron star matter does not 
occur. 

In the present work, we study the properties of polarized neutron matter at finite tem- 
perature using the LOCV technique employing the AVis potential. 



II. LOCV FORMALISM FOR POLARIZED HOT NEUTRON MATTER 

We consider a system of interacting neutrons with A^^i spin up and N2 spin down 
neutrons. The total number density (p) and spin asymmetry parameter (6) are defined as 

p = pi + P2, 

^ ^Vl iV2 

N ^ ' 

5 shows the spin ordering of matter which can have a value in the range of 5 = 0.0 (unpo- 
larized matter) to 5 = LO (fully polarized matter). 

Now, we consider a trial many-body wave function of the form 

^ = F<P, (2) 

where is the uncorrelated ground state wave function (simply the Slater determinant of 
plane waves) of independent neutrons and F = F{1 ■ ■ ■ N) is an appropriate N-body 
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correlation operator which can be replaced by a Jastrow form i.e., 



i>j 



(3) 



in which S is a symmetrizing operator. We consider a cluster expansion of the energy 
functional up to the two-body term, 

1 {i-\H4,} 



For hot neutron matter, the one-body term Ei is 

1=1,2 



El + i?2' 



(4) 



(5) 



Labels 1 and 2 are used instead of spin up and spin down neutrons, respectively, and is 

(6) 



k 



where ni{k,T,pi) is the Fermi-Dirac distribution function, 

ni{k,T,pi 



(7) 



g/3[ei(fc,T,ft)-M,(T,pi)] + 1 

In the above equation P = ^ and pi being the chemical potential, pi is the number density 
and ei is the single particle energy of a neutron with spin projection i. 

In our formalism, the single particle energy, ej, of a neutron with momentum k and spin 
projection i is approximately written in terms of effective mass as {29, 30| 

h^k'' 



ei{k,T,pi 



+ U^{T,pi)- 



(8) 



2m*(p,T) 

In fact, we use a quadratic approximation for single particle potential, incorporated in the 
single particle energy as a momentum independent effective mass. f/i(T, pi) is the momentum 
independent single particle potential. The effective mass, m*, is determined variationally 
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41 



42 



4J. 



The chemical potentials, /ij, at any adopted values of the temperature (T), number density 
(pi) and spin polarization (5), are determined by applying the constraint. 



^ni{k,T,pi) = Ni' 

k 



(9) 



This is an implicit equation which can be solved numerically. 



The two-body energy E2 is 



2A 



(10) 



where 

u{12) = -|i[/(12),[V?2,/(12)]] + /(12)\/(12)/(12), /(12) and V{12) are the two-body 
correlation and potential. For the two-body correlation function, /(12), we consider the 



following form 44 



45|: 



/(12) = 5:/(^)(12)0('=)(12), 



(11) 



k=l 



where, the operators 0^''\l2) are given by 



(12) 



and 5*12 is the tensor operator. After doing some algebra, we find the following equation for 
the two-body energy: 



E2 



^ J^^±^[l_(_l)^+5+ii 



_2_ fh^ 

^vUmJ,^^ 2(25 + 1) 



11 I CO 



dr 



+ 7^({^= - + K - War + 2{Vt - 3Kr) + 2K.}a« (A;^r) 

+ [V^2-3V,2. + V^2.-3Vl2..]c«'(A:/r))(/<i'^)']+ E V) 

+7^({^= + K + K + Kr + (-6A; + U){Vu + - (A; - 1)(V^.. + V^.) 

+ [VV + + (-6A; + 14)Vir][2 + 2Vr,]}a^^^\kfr) 

+ [Vn + 1^2. + Vi2r + V«2..]ci'=)'(A;/r) + [yis2 + V^.2.]rf(f )'(A:/r))/«'^ 

+ ^{V^, + Vlsr - 2{Vn + 1^2. + Vl2.r + Vi2.) " 3(^s2 + ^s2r ) ( A^/r) /f 



(13) 



where a = {J, L, S, Sz} and the coefficient a^^^^, etc., are defined as 



a«7r,p,T)=r^Jz.,5.(r,p,r), 



(14) 



a(,2)^(r, p, T) = r^[/3/j_i,5.(r, p, T) + 7/j+i,5.(r, p, T)], 



a(f)'(r, p, T) = r^[7/j-i,& (r, p, T) + /5/j+i,&(r, p, T)], 



(15) 
(16) 



c^f) ^ (r, p, T) = [r/3/j-i,5. (r, p, T) + udj+^,s. (r, p, T)] , 



'(r, p, T) = r2[e2/j„i,5.(r, p, T) + A2/j+i,5.(r, p, T)], 



2 (r, p, T) = r2[e3/j_i,5.(r, p, T) + A3/j+i,5.(r, p, T)], 



with 



2J + 1' ' 2J + 1' 



2J + 1 



J'(J + 1) j3^2j2 + 3J + 2 

vi = L[L + 1), z/2 = ^ ^ , , , z/3 



2J + 1 



2J+ 1 



^72 



J(J2 + 2J+1) J(j2 + j + 2) 



2J + 1 



^3 



2J+ 1 



J3 + 2 J2 + 2 J + 1 J( J2 + J + 4) 

« — o T , 1 ' S3 — 



2J+ 1 



2J + 1 



and 



_J(J2 + J + 1) j3_^2j2 + 5J + 4 

A2 — „ T , 1 ; A3 



2J+ 1 



2J+ 1 



Ij,sAr,P,T) 



1 



dkidk2ni{ki,T, pi)nj{k2,T, pj)j'j{\k2 - ki\r)- 



27r6p2 

In the above equation Jj(x) is the Bessel's function . 

Now, we minimize the two-body energy Eq. (fT3!) . with respect to the variations 
functions fa^^\ but subject to the normahzation constraint 45 1. 



In the case of spin polarized neutron matter, the function hs^{r) is defined as 

-1/2 



. P . 



] S, = ±l 



; = 



(30) 



where 



riiik^ T, pi)Jo{kr)k'^dk- 



(31) 



From the minimization of the two-body cluster energy, we get a set of coupled and uncoupled 
differential equations which are the same as presented in Ref. j^sj], with the coefficients 
replaced by those indicated in equations f|T4|) - fl22|) . We can obtain correlation functions by 
solving the differential equations and then the two-body energy is computed. 

Finally, we must calculate the total free energy per particle, F, to get the macroscopic 
properties of hot neutron matter. 



F = E-TS, 



(32) 



where S is the entropy per neutron, 

S{p,T) = Zlii^ " ni{k,T,pi)]ln[l - ni{k,T,pi)] + ni{k,T, pi)lnni{k,T, pi)}. (33) 

1=1,2 k 



We introduce the effective masses, m*, as variational parameters [39|, |40|, |41|, |42|, |43|]. In 
fact, we minimize the free energy with respect to the variations in the effective masses and 
then, we obtain the chemical potentials and the effective masses of spin-up and spin-down 
neutrons at the minimum point of the free energy. This minimization is done numerically. 



III. RESULTS 



In Fig. 1, we have plotted the effective mass of spin- up and spin-down neutrons versus 
spin polarization(5) for fixed density p = 0.16 fm~'^ at T = 10 and T = 20 MeV. We see 
that the effective masses of the above components have the same values at 5 = and get 
different values by increasing the spin polarization. It can be seen that the effective mass 
of spin-up neutrons is larger than spin-down neutrons for 6 > and the effective mass of 
spin-down (spin-up) neutrons decreases (increases) by increasing the polarization. In Fig. 
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1, we have also included the results of BPRRV calculations 



31| for comparison. It is seen 



that the behavior of the effective mass versus 6 obtained from our method shows the same 
properties as BPRRV calculations [31 1. 

Our results indicate that the effective masses of spin-up and spin-down neutrons increase 
by increasing the temperature. Whereas, the results of BPRRV jsil and Isayev 46 1 show 
that the effective masses of nucleons in the spin polarized neutron matter and nuclear matter 



are decreasing functions of the temperature. However, in the references [3l|] and 46|], the 
authors determine the momentum independent effective masses through the derivatives of 
the self-consistently calculated single particle potentials at the corresponding Fermi momenta 
while, in our calculations, the effective masses are introduced in the single particle energies 
and considered as variational parameters for the free energy. We note that increasing the 



effective mass by increasing the temperature can be also seen in the resu 



40 



ts of others where 



41 



43|. 



the procedures for finding the effective mass are the same as ours [39|, 

The behavior of the free energy per particle of the polarized hot neutron matter versus 
total number density (p) for different values of the spin polarization (6) at T = 10 and 
T = 20 MeV is shown in Fig. 2. This figure shows that the free energy increases with 
increasing both density and polarization. It is seen that the minimum value of free energy 
occurs for unpolarized case {6 = 0.0) at any density and at relevant considered temperature. 
By comparison, in the two panels of Fig. 2, we see that the free energy decreases by 
increasing temperature. It is also seen that there is no crossing of the free energy curves for 
different polarizations. Conversely, by increasing the density, the difference between the free 
energy of neutron matter at different polarizations becomes more appreciable. According 
to this result, the spontaneous phase transition to a ferromagnetic state in the hot neutron 
matter does not occur. If such a transition existed, a crossing of the energies of different 
polarizations would have been observed at some density, indicating that the ground state of 
the system would be ferromagnetic from that density on. 

In Fig. 3, we have presented the free energy per particle as a function of the quadratic 
spin polarization (5^) for fixed density p = 0.36/m^'^ at T = 10 and T = 20 MeV. It can 
be seen that the free energy per particle increases by increasing the polarization. We see 
that the variation of the free energy of hot neutron matter versus 5^ is nearly linear. This 
indicates that the ground state of hot neutron matter is paramagnetic. In Fig. 3, we have 
compared our results with those of BPRRV calculations 311] ■ It is seen that there is an 
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overall agreement between our results and the results of BPRRV 31|- 

The magnetic susceptibility, Xi which characterizes the response of a system to the mag- 
netic field is given by 

X = -tJ^, (34) 



Js=o 

where /i is the magnetic moment of neutron. In Fig. 4, we have shown the ratio Xf/x 
a function of the total number density at two values of temperature T = 10 and T = 20 
MeV. xf is the magnetic susceptibility for a free Fermi gas. As can be seen from Fig. 
4, even at high densities, the ratio Xf/x increases monotonically and continuously as the 
density increases for any temperature. This shows that there is no magnetic instability in 
hot neutron matter. In Fig. 4, the results of BPRRV calculations 3l|] at T = 20 MeV are 
also presented for comparison. We see that our results and BPRRV results [311 have an 
agreement at low densities. 

In Fig. 5, the difference of the entropy per particle of fully polarized and unpolarized 
cases is plotted as a function of the total number density at T = 10 and T = 20 MeV. Fig. 
5 shows that for all relevant densities, this difference has negative values. According to this 
result, we can conclude that the fully polarized case is more ordered than the unpolarized 



case. In Fig. 5, the results of BPRRV calculations [3l|] are also given for comparison. There 
is an agreement between our results and the BPRRV [3l| results at low densities. 

In Fig. 6, we have plotted the entropy per particle of hot neutron matter versus spin 
polarization for fixed density p = 0.32/m^^ and temperature T = 20MeV. It is seen 
that the entropy decreases as polarization increases with its highest value occurring for 
the unpolarized case. To prevent anomalous behavior of entropy as a function of spin 
polarization, for a given density and temperature, a condition for effective masses can be 



derived 



29|, 



^*(P^^ = ^-^l < 2^/3, (35) 
m*{p,6 = 0.0) ' ^ ^ 

where m*{p,5 = 1.0) and m*{p,6 = 0.0) are the effective masses of the fully polarized and 
unpolarized neutron matter, respectively. From Fig. 1, we have found out that for density 
p = 0.16/m-3, the above ratio at T = 10 MeV is Z'''t7=o^o) = 1-22 and at T = 20MeV is 
"!^^f7^n°n^ = 1-18 which are smaller than the indicated limit. This condition is satisfied for 

mj 2(0=0.0) 
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all other densities explored in this work. Therefor, we can see that the entropy of polarized 
case is always smaller than the entropy of unpolarized case. 

The equation of state of hot polarized neutron matter, P{p, T, S), can be simply obtained 
using 

P(,,r,i) = p=5^(|M (36) 

In Fig. 7, we have presented the pressure of neutron matter as a function of the total number 
density (p) for different polarizations at T = 10 and T = 20 MeV. We see that the equation 

of state becomes stifTer by increasing the polarization. 

IV. SUMMARY AND CONCLUSIONS 

Some thermodynamic properties of the polarized neutron matter at finite temperature 
were reexamined using the lowest order constrained variational (LOCV) method employing 
the ^1^18 nucleon-nucleon potential. Our main goal was to check the occurrence of the 
spontaneous transition to the ferromagnetic state. We found no indication for the occurrence 
the ferromagnetic phase, in agreement with the results of others who used the different many- 
body techniques. Effective mass, free energy per particle, magnetic susceptibility, entropy 
per particle, and the equation of state for the polarized neutron matter at finite temperature 
were calculated and the effect of polarization on these properties were examined. 
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FIG. 1: Our results for the effective mass of spin-up and spin-down neutrons versus spin polarization 
(S) for density p = 0.16 fm.-^ at T = 10 and T = 20 MeV. The results of BPRV calculations [31] 
are also given for comparison. 
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P (fm-3) 




FIG. 2: The free energy per particle of polarized hot neutron matter as a function of the total 
number density (p) for different values of the spin polarization (6) at T = 10 (a) and T = 20 MeV 
(b). 
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FIG. 3: Our results for the free energy as a funetion of the quadratic spin polarization ((5^) at 
r = 10 MeV (full curve) and T = 20 MeV (dashed curve) for p = O.Se/m"^. The results of BPRV 
[31] (dashed curve) are also given for comparison. 
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FIG. 4: The magnetic susceptibility of the hot neutron matter versus total number density (p) 
at T = 10 MeV (full curve) and T = 20 MeV (dashed curve). The results of BPRV [31] (dotted 
curve) are also given for comparison . 




FIG. 5: As Fig. 4 but for the entropy difference of fully polarized and unpolarized cases . 
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FIG. 6: Our results (full curve) for the entropy per particle as a function of the spin polarization 
(S) at T = 20 MeV and p = 0.32/m-^. The results of BPRV [31] (dashed curve) are also given for 
comparison. 
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FIG. 7: The equation of state of the hot neutron matter for different values of the spin polarization 
(6) at r = 10 (a) and T = 20 MeV (b). 



18 



